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We use data collected by the CLEO III detector at CESR to measure the inclusive yields of D s 
mesons as B(T(5S) -> D S X) = (44.7 ± 4.2 ± 9.9)% and B(T(4S) -> D S X) = (18.1 ± 0.5 ± 2.8)%. 

From these measurements we make a model dependent estimate of the ratio of _Bs*'_B a ' ' to the 
total bb quark pair production of (16.0 ± 2.6 ± 5.8)% at the T(5S) energy. 
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PACS numbers: 13.25.Hw, 13.66.Bc 



An enhancement in the total e + e~ annihilation cross- 
section into hadrons was discovered at CESR long ago 
0,0 and its mass measured as 10.865±0.008 GeV. This 
effect was named the T(5S) resonance. Theoretical mod- 
els 0,0,0 predict the different relative decay rates of the 
T(5S) into combinations of B^~B { * ] and B^B^ where 
(*) indicates the possible presence of a B* meson. This 
original ^116 pb _1 of data failed to reveal if B s mesons 
were produced. It is important to check the predictions 
of these and other models; furthermore, e + e~ "B facto- 
ries" could exploit a possible B s yield here as they have 
done for B mesons on the T(4S). 

In this Letter we examine D s yields because in a sim- 
ple spectator model the B s decays into the D s nearly all 
the time. Since the B — > D S X branching ratio has al- 
ready been measured to be (10.5 ± 2.6 ± 2.5)% 0, we 
expect a large difference between the D s yields at the 
T(5S) and the T(4S) that can lead to an estimate of the 
(*) — (*) 

size of the Bs B s component at the T(5S). When we 
discuss the T(5S) here, we mean any production above 
what is expected from continuum production of quarks 
lighter than the b at an e + e~ center-of-mass energy of 
10.865 GeV. The CLEO III detector is equipped to mea- 
sure the momenta and directions of charged particles, 
identify charged hadrons, detect photons, and determine 
with good precision their directions and energies. It has 
been described in detail previously in references and 


In this analysis we use 0.42 fb 1 of integrated luminos- 
ity taken at the T(5S) peak in Feb. 2003. We also use 
6.34 fb -1 of integrated luminosity collected on the T(4S) 
and 2.32 fb -1 of data taken in the continuum 40 MeV in 
center-of-mass energy below the T(4S). These data were 
accumulated between Aug. 2000 and June 2001. The 
detector hardware wasn't changed over the entire time 
period. Efficiencies are carefully monitored and did not 
change measurably between data sets. 

We look for D s candidates through the reconstruc- 
tion of three charged tracks in hadronic events via the 
Dj — > 4>it + decay mode. Here and elsewhere in this 
paper mention of one charge implies the same consid- 
eration for the charge-conjugate mode. Requiring the 
Fox- Wolfram shape parameter i?2 li| to be less than 0.25 
suppresses continuum background events which are less 
isotropic than 6-quark events. 

Pairs of oppositely charged tracks were considered can- 
didate decay products of a <fi if at least one of the tracks 
is identified as a kaon, and if the invariant mass of the 
K + K~ system is within ±10 MeV/c 2 of the nominal 
(f> mass. A third track was combined with the K + K~ 
system to form a D s candidate without using particle 
identification. 

The Ring Imaging Cherenkov (RICH) of the CLEO 



III detector is used for track momenta larger than 
0.62 GeV/c. Information on the angle of the detected 
Cherenkov photons is translated into a Likelihood of a 
given photon being due to a particular particle analyzed 
with a specific mass hypothesis. Contributions from all 
photons associated with a particular track with one mass 
hypothesis are then summed to form an overall Likeli- 
hood denoted as £, for each "z" particle hypothesis. 

To utilize the information on the ionization loss in the 
drift chamber of the CLEO III detector, dE/dx, we cal- 
culate the differences between the expected and the ob- 
served ionization losses divided by the error for the pion 
and kaon hypotheses, called ov and ok ■ 

We use both RICH and dE/dx information in the fol- 
lowing manner: (a) If neither RICH nor dE/dx informa- 
tion is available, then the track is accepted, (b) If dE/dx 
is available and RICH is not, then we insist that kaon 
candidates have PID dE — C 2 - - a 2 K > 0. (c) If RICH 
information is available and dE/dx is not available, then 
we require that PIDrich = —2 log^^) + 2 log(£j<-) > 
for kaons. (d) If both dE/dx and RICH information are 
available, we require that {PIDdE + PI D rich) > for 
kaons. 

To suppress combinatoric backgrounds, we take advan- 
tage of the polarization of the cj> as it is a vector particle 
while the other particles in this decay are spinless. The 
expected distribution from real <f> decays varies as cos 2 dh, 
where 8h is the angle between the D s and the K + mo- 
menta measured in the <j) rest frame while combinatoric 
backgrounds tend to be flat. Thus, we require | cos^l to 
be larger than 0.3. 

For (j)TT + combinations satisfying the previous require- 
ments, we look for D s candidates having a momentum 
less than half of the beam energy. Instead of momentum 
we choose to work with the variable x which is the D s 
momentum divided by the beam energy, to remove differ- 
ences caused by the change of the beam energies between 
continuum data taken just below the T(4S), at the T(4S) 
and at the T(5S). The (fnr invariant mass distributions for 
x < 0.5 are shown in Fig.^ 

We fit the invariant mass of the (pn^ 1 candidates in 10 
different x intervals (from to 0.5) for data taken at the 
T(4S) peak, at the continuum below the T(4S) and at 
the T(5S) peak. 

The invariant mass distribution in each x interval of 
the T(4S) data set was fit to a Gaussian signal shape 
and a linear background. The width of each Gaussian 
was allowed to float. The corresponding distributions at 
the other energies were similarly fit, but with the corre- 
sponding Gaussian widths fixed to those determined at 
the T(4S). The raw D s yields are listed in the second, 
third and fourth columns of Table [I] 

The number of D s candidates is determined by sub- 
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FIG. 1: The invariant mass distributions of the D 3 candidates 
with x < 0.5 from: (a) the T(4S) on-resonance data (b) the 
continuum below the T(4S) resonance data (c) the T(5S) on- 
resonance data. 



tracting the scaled four-flavor (u, d, s and c quark) con- 
tinuum data below the T(4S) from the T(4S) and from 
the T(5S) data. To determine the scale factors, S n s, 
we account for both the ratio of luminosities and the s 
dependence of the continuum cross section using 



E 



nS 



(i) 



where L n s, L cont , E n s and E cont are the collected lu- 
minosities and the center-of-mass energies at the Y(nS) 
and at the continuum below the T(4S). We find: S4S = 
2.712 ± 0.001 ± 0.043 and S 5S = (17.15 ± 0.01 ± 0.24) • 
10~ 2 . 

The second (systematic) error on these scale factors is 
determined by using the number of charged tracks in the 
0.6 < x < 0.8 interval. The lower limit is determined 
by the maximum value tracks from BB events can have, 
including smearing due to the measuring resolution, and 
the upper limit is chosen to eliminate radiative electro- 
magnetic processes. Since the tracks should be produced 
only from continuum events, we suppress beam-gas and 
beam- wall interactions, photon pair and r pair events us- 
ing strict cuts on track multiplicities, event energies and 
event shapes. (Since particle production may be larger 
at the higher T(5S) energy than the continuum below 
the T(4S I ), we apply a small multiplicative correction of 
(0.6±1.1)%, as determined by Monte Carlo simulation to 
the relative track yields.) We find that the scale factors 
using this track counting method are 2.668±0.007, and 
(17.085 ±0.207) • 10 2 , for S 4S and S 5S , respectively, and 
use the difference as the systematic error. 

The total number of hadronic events above four-flavor 
continuum are N^ s) equals (6.43 ± 0.01 ± 0.41) ■ 10 6 



and N$?g s) equals (0.130 ±0.001 ±0.022) ■ 10 6 . The 6.4% 
and 17.5% systematic errors here are due to the 1.6% 
and 1.4% systematic errors on S4S and S$s scale factors 
respectively. 

The branching ratio of T(nS) — > D S X in each i-th x 
interval is given by 



B*(T(nS) -> D S X) 



1 



JV T(nS) 



N*™ s) ■ B(D S -» <Sm) ■ B{<j> -» K+K- 



(2) 



where Ni., C \ are the continuum subtracted on resonance 
D s yields. B(<f> -> K+K~) is taken as 49.1% g. The 
reconstruction efficiency e l is taken to be the same at 
both resonances. This is reasonable because our tracking 
and particle identification efficiencies are carefully moni- 
tored and did not change significantly between data sets. 
Specifically, our Monte Carlo simulations of the D s re- 
construction efficiencies includes time dependent effects 
of dead channels and individual hit efficiencies in both 
the tracking and RICH systems. A comparison of the 
simulations at both energies shows changes in the recon- 
struction efficiency between the two energies of <2%. 

The results are listed in Table ||| We show in Fig. [21 
the x distribution of the inclusive D s yields from T(4S) 
and T(5S) decays, continuum subtracted and efficiency 
corrected. 
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FIG. 2: Branching rate as a function of x from T(5S) decays 
(filled circles) and from T(4S) decays (open circles). 



The total production rate is found by summing the 
partial production rates. The product of the D s produc- 
tion rate at the T(4S) and the B(D S — » 4>ir) is 



B(T(4S) - D S X)-B(D S -> cjm) 



.0±0.2±0.9)-10~ 3 
(3) 
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AeE^I ONT ( 4S ) ONT(5S) Continuum N* r(4S) N^sj e\%) B\ s {%) Bj s (%) 

0.00-0.05 44.4 ±15.7 1.0 ± 3.2 0.0 ± 0.0 44.4 ±15.7 1.0 ± 3.1 28.9 0.11 ± 0.04 0.1 ± 0.4 

0.05-0.10 317.6 ±39.6 13.3 ±8.1 20.7 ± 12.0 261.4 ±51.2 9.7 ± 8.3 23.9 0.8 ± 0.2 1.4 ±1.2 

0.10-0.15 583.6 ± 53.9 30.4 ± 10.4 21.6 ± 15.3 524.9 ±68.1 26.7 ± 10.7 24.7 1.5 ± 0.2 3.8 ±1.5 

0.15-0.20 845.5 ± 59.0 54.4 ± 13.0 41.7 ± 18.5 732.3 ± 77.5 47.2 ± 13.3 25.4 2.1 ±0.2 6.5 ± 1.8 

0.20-0.25 1206.4 ±60.6 57.6 ± 12.7 40.2 ± 18.3 1097.4 ±78.2 50.7 ± 13.0 27.7 2.8 ± 0.2 6.4 ±1.7 

0.25-0.30 2028.6 ±63.8 104.1 ± 14.0 70.3 ± 18.0 1838.0 ± 80.3 92.0 ± 14.3 28.6 4.6 ± 0.2 11.3 ±1.8 

0.30-0.35 2233.7 ±60.7 86.7 ±12.1 57.0 ± 16.2 2079.2 ± 74.9 76.9 ± 12.4 29.4 5.0 ± 0.2 9.2 ±1.5 

0.35-0.40 660.8 ± 37.9 53.8 ± 9.4 75.0 ± 14.5 457.4 ± 54.6 41.0 ±9.7 30.4 1.1 ±0.1 4.7 ±1.1 

0.40-0.45 233.5 ± 25.9 22.6 ± 6.7 73.4 ± 12.9 34.3 ± 43.3 10.1 ± 7.0 31.4 0.1 ±0.1 1.1 ± 0.8 

0.45-0.50 245.8 ±22.2 14.8 ± 5.6 86.0 ± 12.1 12.6 ± 39.5 0.1 ± 6.0 32.4 0.03 ± 0.09 0.0 ± 0.6 



TABLE I: The x dependent D s yields from the T(nS) data, the continuum below the T(4S), the T(nS) continuum subtracted 
data, A^yj^s), the efficiency e 1 , and the partial branching ratios B l nS =T{nS) — > D S X, for ns equal to 4S and 5S. The errors 
are statistical only. 



which is in a good agreement with previous measure- 
ments 0, while at the T(5S) 

B(T(5S) -> D S X)-B(D S -> (/)it) = (19.8±1.9±3.8)-1(T 3 . 

(4) 

Many systematic errors cancel in the ratio of decay rates. 
Thus 



g(T(5S)->j?,X) +0 . 8 
B(T(4S) -> £> S X) U '° 1 ! 



0.3 ' 



(5) 



directly demonstrates, at 5.6 standard deviation signifi- 
cance, a much larger yield of D s at the T(5S) than at 
the T(4S). 

We use B{D S -> </>tt + ) = (4.4 ± 0.5)%, which is the 
weighted average of the (3.6 ± 0.9)% PDG value [6| and 
the recent measured value of (4.8 ± 0.6)% 10], although 
the latter value is at the 90% c.l. upper limit found 
previously [n). We find 



B(T(4S) D a X) = (18.1 ± 0.5 ± 2.8)%, (6) 
and consequently: 

B(B —* D S X) = (9.0 ± 0.3 ± 1.4)% . (7) 
In addition, we find 

B(T(55) -> D S X) = (44.7 ± 4.2 ± 9.9)% . (8) 

From these results, we estimate the size of B' S 'B\ ' 
component at the T(5S) in a model dependent manner. 
Here we start with the knowledge that an equal admix- 
ture of B° and B + mesons decay into the sum of D° 
and D + mesons roughly 100% of the time jfj. Thus we 
expect B s mesons to decay into D s mesons also about 
100% of the time. In what follows we estimate our own 
theoretical corrections to this number. 

We know that the branching fraction B(B — > D S X) — 
(9.0±0.3±1.4)% comes either from the W~ — > cs process, 
shown in Fig. [3J a), or from the b — > c piece if it manages 
to create an ss pair through fragmentation, see Fig.yjb). 




FIG. 3: Dominant decay diagrams for a B meson into D a 
mesons (q is either a u or d quark). 




FIG. 4: Dominant decay diagrams for a B a meson into D a 
mesons. 



Similarly, the production of D s mesons from B s decay 
arises from two dominant processes. Fig. Ufa) shows the 
simple spectator process that is expected to produce D s 
mesons nearly all the time; here the primary b — > c tran- 
sition has the charm quark pairing with the spectator 
anti-strange quark. Fig. Efb) shows the subset of pro- 
cess (a) where W~ — » cs and these two quarks form a 
color singlet pair. The chances of this occurring should 
be similar to the chance of getting an upper- vertex D s in 
B decay (Fig.^a)), i. e. a D s along with a D. 

We can use data to help estimate the size of these 
processes. First let us consider the diagram shown in 
Fig. E[a). The nearly 100% probability that this pro- 
cess will produce D s mesons is reduced if the cs pair 
fragments into a kaon plus a D instead of a D s by pro- 
ducing an additional uu or dd pair. We don't actually 
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know the size of this fragmentation, though it's clear 
that producing a light quark-antiquark pair (dd or uu) 
is easier than ss. We estimate that the reduction in 
D s yield due to this fragmentation is a (—15 ± 10)% ef- 
fect. Next we estimate the size of the process depicted 
in Fig. IHb) . The B —* DD S modes have branching frac- 
tions that sum to about 5%. There are some additional 
decays due to B — > D**D S and B — » DD^*] decays, 
that also contribute D s mesons. We add these and es- 
timate an extra (7 ± 3)% of D s mesons in B s decays 
produced by diagram Fig.QJb). Taking into account all 
these contributions, we derive a model dependent esti- 
mate of (100 + 7 - 15)% = 92%. Therefore, we use 
B(B S -> D S X) = (92 ±11)%. 

We can estimate now the fraction of the T(5S) that 
(*) — (*) 

decays into B s B s , which we denote as f s . The D s 
yields at the T(5S I ) come from two sources, B and B s 
mesons. The equation linking them is 



and B(T(4S) -> D S X) = (18.1 ± 0.5 ± 2.8)%. The ratio 



B(T(55) -> D S X)B{D S -> <M + )/2 = 
f s -B{B s ^D s X)B{D a ^<t>iT + ) 
(1 - /.) 



B(T(45) -> D S X)B(D S -> cfm + ) 



where the product branching fractions B(T(5S) 



(9) 



£(£ s 



6tt+) and 6(T(4S*) 



D S X) 



D S X) 

B(D S 07T+) are given by equations 0] and OH respec- 
tively. Therefore, at the T(5S) energy, we obtain the 

(*) — (*) • — 
Bs B s ratio to the total bb quark pair production 

above four-flavor (it, d, s and c quarks) continuum of 

/, = B(T{5S) -» B^B^) = (16.0 ± 2.6 ± 5.8)% (10) 

The systematic errors in this analysis are dominated 
by the 1.6% relative error on S4S and 1.4% on S$s scale 
factors which contribute large components (6.4% and 
17.5%) to the error on the number of hadronic events 
above continuum at the T(4S) and T(5S). There is also 
a contribution from the 11.3% error on the B s — > D S X 
branching fraction estimate and a contribution from the 
11% error on the absolute D s — > 07r branching fraction. 
An additional component comes from a 6.4% error on the 
D s detection efficiency, which includes a 2% error on the 
tracking efficiency and a 2% error on the particle identi- 
fication, both per track. We also have 5% error on the 
yields due to the fitting method. The total systematic 
error is obtained by summing all entries in quadrature. 

In conclusion, we have measured the inclusive yields of 
D s mesons as B(T(5S) -> D S X) = (44.7 ± 4.2 ± 9.9)% 



B(T(5S) -> D S X) 



2.4 ±0.3: 



(11) 



B(T(4S) -> D S X) —-—-0.S1 
provides the first statistically significant evidence (5.6<r) 
of substantial production of B s mesons at the T(5S) res- 
onance. Using a model dependent estimate of B(B S — > 

f*l — (*) — 
D S X), we find that the Bs B s ratio to the total bb 

quark pair production above the four-flavor (u, d, s and 

c) continuum at the T(5S) energy is 



f s = B(T(55) -» B^B { * y ) = (16.0 ±2.6 ±5.8)%. (12) 



Several phcnomcnological models predict the decay 

rates of the T(5S) into combinations of B^B * and 

Bs B s , though here we are only concerned with the 
relative EL fraction f s . The unitarized quark model es- 
timates [J and the predictions of Martin and Ng Q are 
about 30%, both somewhat larger than our measurement. 
Byers and Eichten [f| present two models both giving 
f 8 < 20% in good agreement with our finding. 

We gratefully acknowledge the effort of the CESR staff 
in providing us with excellent luminosity and running 
conditions. This work was supported by the National 
Science Foundation and the U.S. Department of Energy. 



[1] 
[2] 
[3] 



[4] 
[5] 

[6] 
[7] 



[8] 



[9] 

[10] 
[11] 



D. Besson et al. (CLEO Collaboration), Phys. Rev. Lett. 
54, 381 (1985). 

D. M. Lovelock et al. (CUSB Collaboration), Phys. Rev. 
Lett. 54, 377 (1985). 

S. Ono et al, Phys. Rev. Lett. 55, 2938 (1985); S. Ono, A. 
I. Sanda, N. A. Tornqvist, Phys. Rev. D 34, 186 (1986); 
N. A. Tornqvist, ibid 53, 878 (1984). 

A. D. Martin and C.-K. Ng, Z. Phys. C 40, 139 (1988). 
N. Beyers and E. Eichten, Nucl. Phys. B (Proc. Suppl.) 
16, 281 (1990); N. Beyers [hep-ph /9412292 (1994). 

S. Eidelman et al, Phys. Lett.. B592, 1 (2004). 
D. Peterson et al, Nucl. Instrum. and Meth. A478, 142 
(2002); Y. Kubota et al. (CLEO), Nucl. Instrum. and 
Meth. A320, 66 (1992). 

M. Artuso et al, " The CLEO RICH Detector," 
arXiv:physics/0506132 and M. Artuso et al, Nucl. In- 
strum. and Meth. A502, 91 (2003). 

G. Fox and S. Wolfram, Phys. Rev. Lett.. 41, 1581 
(1978). 

B. Aubert et al et al. Phys. Rev. D 71, 091104(R) (2005). 
F. Muheim and S. Stone, Phys. Rev. D 49, 3767 (1994). 



